Sympathetic cooling route to Bose-Einstein condensate and Fermi-liquid mixtures 
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We discuss a sympathetic cooling strategy that can successfully mitigate fermion-hole heating in a 
dilute atomic Fermi-Bose mixture and access the temperature regime in which the fermions behave 
as a Fermi-liquid. We introduce an energy-based formalism to describe the temperature dynamics 
with which we study a specific and promising mixture composed of ^Li and *''Rb. Analyzing the 
harmonically trapped mixture, we find that the favourable features of this mixture are further 
enhanced by using different trapping frequencies for the two species. 
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Introduction. Cold atom mixtures of fermions and Bose- 
Einstein condensates (BEC) suggest a new venue for 
the study of quantum liquid mixtures that were previ- 
ously available only in the form of condensed '^He-^He 
mixtures Ql ■ While we understand the helium phase sep- 
aration to be caused by mediated interactions |3| , a first- 
principle description of this and other helium quantities is 
complicated by strong interaction effects 0. These also 
reduce the boson mediated interactions that can Cooper- 
pair the '^He fermions, giving a critical temperature Tc 
currently unaccessible. Mediated interactions and their 
prominence near a quantum phase transition (quantum 
criticality, see 0) have also been identified as the prob- 
able cause of the observed non-Fermi-liquid behavior in 
high Tc superconductors and, perhaps, of high su- 
perconductivity itself. In these systems, however, a clear 
understanding is again impeded by the ambiguity of ex- 
perimental data. 

The above issues can be studied in a much cleaner 
and more accessible environment when cold atom physi- 
cists succeed in cooling fermion-boson mixtures into the 
Fermi-liquid regime |2| in which the fermions behave in 
a universal manner that follows from expanding the en- 
ergy to second order in quasi-particle occupation num- 
bers. This implies a heat capacity that varies linearly 
with temperature, a feature that often serves as experi- 
mental test of Fermi-liquid behavior. In a dilute fermion 
gas, it sets in below ten percent of the Fermi tempera- 
ture T < O.lTp, whereas T < 0.02Tp is required in the 
■^He liquid. For cold atom quantum liquid studies, which 
can also access fermion-boson mixtures in a parameter 
regime that is fundamentally different from that of the 
helium mixtures , crossing into the Fermi- liquid regime 
is as crucial of a step as cooling below Tc is for the study 
of superfluidity. So far, experiments that use bosons to 
cool fermion atoms (sympathetic cooling) have not en- 
tered this regime because the low temperature drop 
in the BEC heat capacity impairs its ability to refriger- 
ate Isl , and because the unavoidable loss of fermions leads 



to significant heating The BEC introduces further 
losses in the Fermi system (three-body recombinations) 
thereby increasing the Fermi-hole heating rate. 

In this Rapid, we show that with a large boson-to- 
fermion mass ratio and with the proper trapping param- 
eters, sympathetic cooling can access the Fermi-liquid 
temperature regime. Heavy bosons require many colli- 
sions from the lighter fermions before they heat signifi- 
cantly. The corresponding increase in BEC heat capacity 
is crucial; low three-body recombination rates are help- 
ful. The need to strike a delicate balance between keeping 
the fermion boson overlap sufRciently low to reduce hole- 
heating, but sufficiently high to ensure efficient fermion- 
BEC heat exchange can be met by adjusting the fermion 
and boson trapping frequencies in a bichromatic trapping 
scheme. We illustrate these points for the promising can- 
didate of a ^Li-®^Rb mixture, which has been recently 
experimentally investigated 

Loss-induced Temperature Dynamics. - The cold atom 
many-body states of interest are metastable (the true 
ground state being a condensed solid) and trapped atom 
gases have a finite lifetime. Atoms leave the trap as a 
consequence of collisions with room temperature back- 
ground atoms, of two-body spin relaxation processes (in 
magnetic traps), and of three-body recombination. In- 
stead of describing the many-body response to particle 
loss, we keep track of the energy-balance and define the 
effective temperature of the system T as the temperature 
at which an equilibrium system with the same number of 
particles has the same energy. 

The above collisions reduce the number of particles of 
type j, Nj, as well as internal energy E, as each loss 
removes the energy of the lost particles. If the collision 
products, on their way out of the trap, collide with the 
remaining atoms, the process also adds an average ki- 
netic energy e^°". We equate the total energy balance 
rate, dE/dt — J2j £'j°''^dNj/dt to the variation of the to- 
tal equilibrium energy of a temperature-relaxed system, 
dE{Nj,T)/dt = {dE/dT)dT/dt + J2j{dE/dNj)dNj/dt 
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FIG. 1: Hole heating in homogeneous Fermi-Bose mixtures. 
The temperature, T, scaled by the Fermi temperature, Tf, is 
depicted for four Fermi-Bose mixtures consisting of fermionic 
®Li or ''"K atoms with '^'^Na or *^Rb as bosonic partners, 
starting at T/Tp = 10~^. For comparison, the hole heating 
curve for pure fermions, assuming the same loss-rate, is also 
plotted (lighter continuous curve). 



and we identify the usual thermodynamic derivatives, 
{dE/dT) ~ C, where C denotes the heat capacity at con- 
stant number of particles, and (dE/dNj) = ^j, where 
fij represents the chemical potential of the j-particles. 
Equating the energy rates and solving for the tempera- 
ture derivative, we obtain the central equation of tem- 
perature dynamics 




coll dNj 
dt 



(1) 



where we have defined F as the free energy that occurs in 
Fermi-liquid theory: F = E — "Y^- j-ijNj (here dependent 
on Nj rather than fj,j). 

Homogeneous ideal gas Fermi-Bose mixtures. - The ho- 
mogeneous mixture of a single component degenerate 
fermion gas and a BEC allows a transparent description 
of hole heating and its mitigation by sympathetic cooling. 
Apart from border effects, the homogeneous mixture can 
be realized by using bluc-detuncd laser light sheets [l^ . 
Atoms in the mixture undergo various three-body loss 
processes. For instance, a fermion can rccombinc with a 
boson into a highly energetic molecule while ejecting a 
second boson, corresponding to 



(2) 



where np and ub represent the fermion and boson densi- 
ties respectively, and a denotes the three-body loss rate 
coefficient. Processes involving three-boson collisions do 
not lead to significant heating and those involving two or 
three indistinguishable fermions are Pauli-inhibited. If, 
as in Eq. Q, the rate is energy-insensitive, the loss of one 
out of Np fermions removes, on average, a kinetic energy 



^Ep, where Ep is the Fermi energy. At the same temper- 
ature the remaining degenerate system of A^f — 1 fermions 
would have an energy Ep less than the initial Fermi sys- 
tem of Np fermions. The net result is an effective energy 
increase of ^Ep per fermion lost (assuming, optimisti- 
cally, e^°" = 0). We apply Eq. (QJ to the homogeneous 
mixture of fermions of mass mp, and an ideal BEC of 
A^B bosons of mass mn. The total heat capacity of an 
an ideal or near-deal gas mixture is the sum of the BEC 
and fermion heat capacities, C = ksN p^Tr'^ /2){T/Tf)[1 + 
iA5/8Tr^/'^)C5/2{lKmB/mp)y^,/T/T;i where ks de- 
notes the Boltzmann constant, Tp is the Fermi- 
temperature, Tp = Ep/ks and C5/2(l) = 1-342. The 
time scale of the dynamics is the lifetime r of the fermion 
system, r~ 
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The term oc \JT jTp in the denominator stems from the 
BEC heat capacity and represents its capability to ab- 
sorb heat released from fermion-hole heating. Note the 
sensitive dependence on mass-ratio illustrated by Fig. 
1, where the temperature of mixtures that start out at 
T = O.OlTp is shown for four possible Fermi-Bose mix- 
tures with stable alkalis (^Li or for the fermionic 
species, ■^'^Na or ^''Rb for the bosonic species). For the 
sake of comparison, we recover the fermion-hole heating 
rate for a pure fermion system obtained in Ref. |lCll | by 
setting rriB ^ 0, corresponding to a BEC-refrigerator of 
vanishing heat capacity (shown by the thin full line). 
Trapped interacting mixtures. - In addition to the mass- 
ratio dependence, the spatial density profiles in the re- 
alistic case of trapped gases affect the temperature dy- 
namics. The fcrmion-boson overlap can be reduced by 
a repulsive inter-species interaction or by the occurrence 
of a true phase separation. In the non-phase separated 
systems, the mutual overlap can be controlled by vary- 
ing the ratio of the fermion to boson trapping frequencies 
(using, for instance, a bichromatic trap |9|). A decrease 
in overlap reduces fermion loss but also reduces the ef- 
ficiency of sympathetic cooling. Which effect wins out, 
depends sensitively on the density profiles and, hence, on 
the interaction parameters. 

To determine the trap densities in the case of a ^Li- 
^''Rb mixture, the most promising candidate shown in 
Fig. 1, we have evaluated the interspecies elastic scat- 
tering lengths for various Li and Rb isotopes. Although 
spectroscopic data for LiRb potential curves are lacking, 
the excellent agreement between ab initio ^3 and our 
parametrized potentials (see Table at large separation 
attest to the accuracy of the curves in that region. How- 
ever, there is a discrepancy between both sets of curves 
near the equilibrium distance. Compared to |l3l | our po- 
tential curves are deeper by about 3% for the singlet and 
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and nF(r) in a Thomas- Fermi approximation are deter- 
mined by (see Molmcr reference in |^) 



TABLE I: Singlet and triplet elastic scattering lenghts (shifted 
and unshifted) for the mixtures of Li and Rb isotopes, in Bohr 
radii. The values of s are found so that ar = ±17 ao, and 
the same scaling was assumed for tlie singlet. The scattering 
lengths were obtained from potential curves with ab initio 
data from [l5j |. joined smoothly to an exponential wall of the 
form ce~^^ at short separations R, and to the long-range form 
-Cfi/R^ -Cs/R^ -Cw/R'-^TAR°'e-^"- at i?=13.5 ao. Here, 
+ stands for the singlet X^Ti^ and triplet a^E"*" molecular 
states, respectively. The values for Cs = 2545 a.u. from [l^ . 
and Cs = 2.34 x 10^ a.u. and Cio = 2.61 x lO'' a.u. from 
fr^ . were used. The parameters of the exchange energy are 
a = 4.9417 a.u. and (3 = 1.1836 a.u. [H, while the constant 
A — 0.0058 a.u. was found by fitting the ab initio data. 



7% for the triplet channel. By scaling the whole potential 
curves so that Vs{R) = sV{R) with s < 1, we can explore 
the effect of having shallower potential curves; the values 
of the scattering lengths change significantly as s varies 
accordingly. In the singlet case, one bound level dis- 
appears, and the scattering length varies between ±oo, 
while the triplet scattering length becomes positive. Re- 
cently, Zimmermann and co-workers [lH determined the 
magnitude of the triplet scattering length of ^Li-^^Rb 
to be IotI = 17^6 0,0- In Table U we present values 
for the singlet and triplet scattering lengths for the un- 
shifted (s = 1) and shifted potentials (so that ax agrees 
with ±17 ao). These values as well as those calculated in 
IQl are tentative and future measurements are required 
to specify the potentials more accurately. Such experi- 
mental feedback can be obtained from either one of the 
isotopes so that we list the results for all combinations, 
regardless of their fermionic or bosonic nature. 



We now study the temperature of a ^ Li - ^^Rb mix- 
ture trapped in an idealized trap: both species are con- 
tained by cylindrical trapping potentials Vp (r) and Vb (r) 
that are of the hard wall type in the transverse direction 
with the same radius Rt, where Rt is much larger than 
the BEC healing length and the average fermion-fcrmion 
distance. In the longitudinal direction, the fermions and 
bosons experience harmonic trapping potentials with dif- 
ferent trapping frequencies Uf and tUh respectively. Apart 
from its simplicity, such a situation reproduces the rele- 
vant heating features of a generic three-dimensional cold 
atom mixture with anisotropic confinement. From Table 
I, we find that some hyperfine states of the ^Li and Rb 
isotopes that can be trapped magnetically give positive 
valued inter-species scattering lengths, thereby favoring 
spatial separation in the trap. The density profiles ubIv) 



?^B(r) = 



nF(r) 



MB - VB(r) - AFBnF(r) 
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-{AiF-^F(r)-AFBnB(r)]}3/2 (4) 



where Ab and Afb denote the boson-boson and fcrmion- 
boson interaction strengths proportional to the corre- 
sponding s-wave scattering lengths, as and apB'- Ab = 
4:Trh?aB/mB and Afb = 27rfi,^aFB (l/wp ± 1/tob). In Eq. 
Q), we tacitly assume temperature independent density 
profiles, a good approximation in the regime of interest. 
Solving Eq. numerically, we find that a decrease in 
confinement strength for the bosons generally increases 
the overlap with the Fermi species. Controlling the rel- 
ative trapping frequencies allows for a variation of the 
fcrmion-boson overlap and, hence, the rate of Fermi-hole 
heating and the efficiency of sympathetic cooling. The ef- 
fect of inter-species interactions is marginal on the Fermi 
gas but is very pronounced for the BEC, and the mutual 
overlap is very sensitive to the interaction parameters. 
The temperature trajectories shown in Fig. 2 are calcu- 
lated for ojf = lO'^s"^, while varying the Bose trap fre- 
quency as shown in the caption. We choose the chemical 
potentials in the initial state to be 1 iJ,K for the Fermi 
gas and 100 uK for the Bose gas. The scattering lengths 
arc chosen as aB=5.8 nm for Rb and aFB=-0.90 nm for 
the Li-Rb interactions. 

We consider the fermions to be subject to loss by back- 
ground scattering, ri-p{x) = —'yn-p{x), as well as by the 
3-body recombination of Eq. Q. While these loss-rates 
arc unknown, we have chosen values that corresponded 
to 7 = 10 Hz and t^^ = ari^g = 4s~^ at the peak 
densities of case LOf/oj^ = 1 of Fig. 2. The behavior 
shown in Fig. 2 is quite robust with respect to various 
choices of the relevant parameters. In addition, we as- 
sumed that a continued evaporative cooling of the bosons 
drains away energy from the system at a realistic rate. In 
the calculations with different fermion and boson trap- 
ping frequencies, we took the evaporative cooling rates 
to be equal and accounted for the low overlap drop in 
evaporative cooling efficiency by introducing a figure of 
merit as described in Ref. [l9|. This allows to intro- 
duce an effective evaporative cooling rate —QeS propor- 
tional to the overlap of the two clouds. The rate of free 
energy change is then identified as sum of three terms: 
dF/dt = [dF/dt)^^ + {dF/dt)^^^^^ - Qcft, where the free 
energy rates have to be calculated from the resulting den- 
sity profiles, 

(dF\ {6n^)y'h^a f ,s . , .,5/3 . , ,,2 
[-dt),,= 5mF ["bW] ' 

(5) 

whereas as in [lOj , the background scattering free energy 
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FIG. 2: Fermi hole heating of harmonically trapped ^Li-®''Rb 
mixtures. The dependence of the degeneracy parameter T/Tp 
(starting from an initial value of T/Tp = 10~^) on time is 
depicted for three cases of ojf/ojb = 1 (dot-dashed curve), 
uJi/uih = 2 (dashed), uof/ujh ~ 10 (continuous). The dotted 
curve represents the intrinsic Fermi hole heating for the same 
loss-rate of the fermions but in the absence of the ^'^Rb atoms. 



increase rate is given by 



dF 



Bkend 



^Q^^f'h^l / ,3 r , ,15/3 

-^^-^jdrMv)]^ . (6) 



Finally, the temperature trajectory follows from Eq. l|T|l. 
in which we use Eqs. ^ and © 



dt 



{dF/dt) 



Bksnd 



Cf(T) + Cb(T) 



(7) 



where Cp and Cb denote the fermion and boson heat 
capacities. The resulting temperature dynamics is de- 
picted in Fig. 2 for three different trapping frequency 
ratios. These results show that the heating rate is miti- 
gated best by using a larger wt/wb-ratio, at least until the 
spatial overlap is decreased due to an excessive spreading 
of the Bose cloud. 

Conclusion. - Our study of loss-induced heating, based 
on an effective equilibration model, reveals that a cold 
atom Fermi-Bose mixture of high boson to fermion mass 
ratio such as ^Li-^^Rb can be maintained in the Fermi 
liquid temperature regime. This result, a consequence of 
the increased BEC heat capacity due to the larger trap- 
ping frequency ratio between fermions and bosons, could 
open up a new avenue for cold atom studies. We also 
discuss a specific cooling strategy which could be soon 
implemented in the experimental investigations going on 
for the specific ^Li-^^Rb mixture, for which anomalous 
heating has been observed 0| . 
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